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mSTRACT 


Flexural response of piles to exfemally applied loads 
is one of the oonplox soil-structure interaction prohlcns. For a 
realistic analysis of the pile prohlem it is necessary to consider 
the proper houndary conditions, axial force variation, variation of 
moment of inertia of the pile and load-deflection relationship of the 
pile. It has boon concluded from the literature review .(Ch-apfOr' 2’) 
that the analysis of an individual pile can be cariiod out at working 
loads using the concept of soil modulus whereas at larger loads the 
plastic resistance developed in the top region has to be taken into 
account. 

In Chapter 3» closed form solutions are obtained for 
axially and laterally loaded tapered piles in layered soils for 
tlifforcnt boundary conditions. In Chapter 4 closed form solutions 
are obtained for individual pile in olasto-plastic soil. 

In Chapter 5> after roviev/ing the existing litcra,ture 

on the probabilistic approaches to soil mechanics problems, the need 

for considering the uncertainties involved has boon stressed. Also ' 

of 

an initia,tion has boon made to bring about the influence /factor of 
safety, nxmibGr of samples end variance' on the probability of 
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HTTRODUCTIOIT 

In foundation engineering the studjr of flexural hehaiTleBix 
of piles forms one of the ooniplex soil-straotare interaction 
prohlems, for a realistic analysis of the flexural behaviour •of' 
a pile, 'it is necessary to know the soil pressure distribution at 
the interface between the pile and the soil, whioh intum is 
governed hy a host of parameters. Hence if one has to take all 
the paxameters into account the analysis becomes extremely complex. 
However with the experience and the data available from the field 
tests and laboratory exp)erinients, some simplifying assumptions can. 
he made to arrive at simple analytical solutions for many oases 
corresponding to field situations. 

On the basis cf study of the v/oaic done by various 
investigaters in this field one oan have a broad idea about the 
qualitative nature of contact pressure at the interface. Ihe 
systematic studies carried out hy Matlock and RLppexger ( 19 57 f 1957a) 
Peck and Pavisson (1958, 1962), Broms (1964a, 1964b) and Ke’risal 
(1965) show that the flexural analysis of piles can be made based, 
on the ooneept of asil modulus at working loads. According to this 
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concept, the soil pressuro, p, ia assumed to he directly 
proportional to the deflection at that point. The Turortrs of 
Biot (1957) Vesio'(l 96 l) deal with the validity of this 
assumption in the analysis of heanis on elastic foundation. 

From the analysis of extensive field and lahoratory data, 

Broms (1964 19^4 h) has shorn, that the concept is quite 

valid for analysing the laterally loaded piles at working 
loads. However, if the soil in the top region is quite loose 
and does not talce much load, then there will he large deflection 
or rotation duo to the applied loads, and the problem becomes 
much more complicated due to non-linear characteristics •f the 
soil and this needs modification in applying the concept of 
soil modulus. 

From the available literature, it is seen that the 
flexural behaviour of the pile, in general, is governed by the 
soil properties in the top region upto a depth of 10 to 15 pile 
diameters from ground level. As the properties of the top 
region may be different from those of the lovex region it is 
very important to obtain solution for flexural response of pile 
in a layered soil. This, also brings out the drawbacks of the 
existing literature, where in either constant or linearly 
varying soil modulus is assume 1 , Recent studies of Davisson 
and Gill (1963) j Heddy and Yalsangkar (1968, 1971 )» Basudhar 
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(1971) <3.eal wiidi tlie, prolalem. of flosural analsrsis of unifom 
oross-seotion pile, at v/orking loads, in layered soil. 

Relatively fen solutions are available for tbe flemral analysis 
of pilos whoso moment of inertia varies vrith depth because tf 
mathematioal complexities. Mostof the solutions use numerical 
techniques for the analysis and with this in view in Chapter 3 
generali25ed solutions are presented for laterally and axially 
loaded tapered piles in layered soil, numerical results are 
presented exhaustively, bringing out qualitatively the influence 
of soil modulus variation, layering and aodal force variation. 

Engineers, in addition to uaider standing the behaviour at 
working loads, are also interested in the ultimate carrying 
capacity of the soil -structure system. In case if laterally 
loaded piles, to determine the marimum load carrying oapacity,it 
is assumed that, the soil fails all along the length of pile or 
that the section of the pile -fails, depending upon the length and 
boundary conditions of the pile at top (Eroms, 1964 a, 1964 b) . 
However many field studies indicate that in reality the situatif»n 
lies somov/hat intermediate, in the sense, only some part of the 
soil at the top may go into plas-tic oondition, but the depth of 
plastdo zono may not he sufficient to induce exsessivo bending 
moments, sd that pile failure takes place, Si> for a realistio 
analysis both elastic and plastic properties have to be oonsidered, 
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as "both states occirc along the embedded length of the pile. 
Studies of Wagner (1955)> Matlook, Eipperger and Pitzgibhon 
( 1956 ), Matlock and Eipperger (l957a)> Mori (19^4)/ Eeddy and 
Yalsan^car (l970), M,B. Koy (l970) and Basudhar (1971 ) snppart 
this view, Eeoently, Domenico Lalli (1971) li-as presented 
approximate solutions for -ftie behaviour of piles of variable 
eross-seotional properties aid head angular restraint, in 
hi-linear soils under both lateral and axial loads ly large 

IN 

deflection Matrix method. In the above case the results are 
presented for a pile of particular dimensions and varying 
lateral loads^^ axial loads ard head angular restraint in dipign*" 
sional forms. 

In Chapter 4j analytioal solutions in closed ferm are 
presaited for a tapered pile in eiasto -plastic soil with geno^ 
ralized loading . As the solutions are obtained jp. closed form 
these can be used for checking the oomrergenoe of numerical 
methods used for solving such problems. Also the results oam 
be directly used in field situations wherein most of the assump- 
tions- made in the analysis are satisfied. Even in oases wher^n. 
the axial force variation and moment of inertia variation may he 
different than Iho one assumed in the analysis, the solutiens 
can be used to get fairly good estimate of flexural resptase* 
Mpngdeal data has been obtained oon^rehensively, to bring about 
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the influence of various parameters qualitatively. Results arc 
presented in non-dimmsional form so that they can he useful for 
design purposes. 

It is well Idown that the soil properties vary from point 
to point even in an a/pparantly uniform soil, duo to the random 
nature of the soil. .Is such in recent times many authors have 
considered the uncertainties in the soil properties in the analysis 
of soil mechanics and foundation engineering prohlema. 

In oha-i^ter 5 sji initiation has hcon made for taking into 
account the uncertainties involved hy adopting the probabilistic 
approach in the ultima,te load, ana-lysis of latora-lly loaded piles. 

The oha-ptcr also contains a review of a.11 the available literature 
concerning probabilistic methods as applied to soil mechanics and 
foundation engineering problems. 



GHijPTER 2 


BEVim OP H'lERY/rURE 


2.1 GENERii 

Tills Chapter deals tvith a rovieTf of the experimental and 
analytical investigations pertaining to flemral hohaviour of piles 
under the action of latera-1 or lateral and axial loads. In the 
review of the experimental investigations^ attention is given mainly 
to the recent studies on instrumented piles, which have added to the 
clear understanding of this oom-plex soil- structure interaction 
problem. Also detailed review is ava-ilable in the works of Davisson 
(i960), Talsangkar (1969) and Basudhar (1970* At the end of the 
Chapter, generalized conclusions are drawn, which form the basis for 
the assumptions made in the sub&&q.uent chapters on the analysis of 
piles. 

2 .2 K:{pmiiEiT;ji htvestigatiois 

For the first time the contact pressure distribution at 
the interface vfas measured by Stobie (1930), while conducting 
eicpariments on behaviour of poles. The study brought out the 
fact, that the passive resistance offcoxid by the pile is conside- 
rably more than the two-dimensional passive x)rGssuros obtained 
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ty Ranlcme’s theory. Sut the measareraen.ts .were-'-very crude. 

Experimental iriTCstigations were conducted by Peagin 
(1957) , Mataa (1939), Silts, Graves and Driscoll (l94S) on 
uninstramented piles. The results of these studies bring out 
the significant influence of repeated loading, batter of the 
pile and stiffness of the pile as affecting the flesural behaviour. 

For the first time, systematic studies to measure the 
soil pressure by using strain gauges v/ere conducted by Loos and 
Breth (l949.) on model steel vortical and batter piles. The soil 
modulus was assumed to vary linearly for the analysis of test 
results and the correlation between the predicted values and 
observed values was seen to be good. The results also revealed 
that, increasing the length of the pile beyond a certain liraLt had 
practically no effect on resisting lateral loads. 

MoCammon and Asoherman (1953) presented the results of the 
tests oonduotod on long hollow piles of length of about 170 feet, 
with an overhang of 80 feet embedded in plastic clay. The results 
indicate that the ola.y acted as an elastic medium and the point of 
maxLinum moment oocured very near to the ground surface for both free 
and restrained piles, 

Gleser (1953)’ has shown frcm tho field tests, that the 
presence of vibrating operations near the piles increases the 


deflections 
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Barter (1953)» while discussing the test results of Gleser, 
reported the results of the tests carried out on monotuhe piles with 
or without acdal loads. The test piles were embedded 80 to 90 feet 
in organic olay. From these test results it was concluded that, the 
deflections get considerably magnified in the presence of arial lo 
in addition to lateral loads. 

Maspn and Bishop (l954) presented the results of the tests 
conducted on a full scale pile embedded in earth fill and subjected 
to lateral loads. The pile was held in vertical position and the sand 
around it was well compacted. Mathematical analysis was presented 
based on soil modulus concept and the results of theoretical analysis 
and observed values were reasonably in good agreement. Mason (1956) 
by continuing the studies at -ftie same site with repetitive loading on 
the piles, reported that 1he repetition of loading increases the 
deflection by about 20 per cent* 

McNulty ( 1956 ) carried out tests on Raymond Standard piles 
and on wooden piles. Three unreinforced concrete Raymond piles were 
embedded in a soil deposit upto 14, 22 and 26 feet. The soil at the 
site consisted of 3 feet of sandy olay at the top below which there 
was fine silty sand. The v;ater table was at the ground surface. It 
was interesting to note that in the absence of aodal load unreinforoed 
Raymond pile failed at relatively low loads as practically there was 
no flemral resistance offered. However, when the tests v/ere conducted 
with axial loads, the piles shovired greater lateral resistance as the 
compressive stresses from the axial load kept the tensile stresses 
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caused "by the lateral load to a mirdimim. 

For the first time most exhaustive and systematic eaqperimentad- 
investigations, hoth in laboratory and field were conducted "by Matlock 
and Eipperger (1956, 1957 a, 1957 h). 

liatlock, Ripperger and Fitzgibhon (1956) conducted a series of 
tests on a, steel pile with 12.75 inches outside diameter and 0.5 inch 
wall thickness. The test pile was driven to a depth of 42 ft. in a 
fresh Water lake deposit. The soil reactions along the embedded length 
of the pile were measured by strain gauges. The test site consisted of 
predominantly ol^, with intermediate sand layers. At first instance 

, ' I 

two senes of tests were made, one of the static load and the other of 
cyclio load. 

The results of the static load tests indicate that the soil 
resistance was plastic for large lateral loads for upper 4 to 6 feet of 
soil. It was also observed that, for cyclic load test, the soil 
resistance in the upper 14 feet of soil was plastic and the soil 
resistance was nearly elastic below 14 feet. This considerable 
decrease in the soil resistance in the tipper region compared to the 
static case can be attributed to the considerable remoulding of the soil 
due to cyclic loading. From the analysis of the test results it was 
also seen that the soil modulus variation was linear in the upper 
region and constant in the lower region. Matlock and Ripperger 
(1957 S') contirrued the tests at the same site by varying the soil 
condition at the top 10 feet of the soil deposit. This test brought 
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out the ■beneficial effect of plsuaing sand or gravel in the top region 
around the pile in resisting lateral loads, 

Matlock and Eipperger ( 1957 '^) perfQiBied studies on model 
aluminim piles embedded in gelatin, Photoelastio studies were 
performed on the laterally loaded piles. It was seen from these 
studies that, major part of the load was dissipated within a distance 
of 5 pile diameters from the face of the pile. It wasalco observed 
that at higher loads plastic resistance developed within a distance of 
1,5 pile diameter, from the face of the pile. It was also concluded 
from these tests that, for a given lateral load the deflections of the 
pile were independent of pile diameter. 

Davisson (i960) reports the lateral load tests carried- out by 
California Division of Hi^ways in 1958, on reinforced concrete piers 
oast in bored holes. The soil consisted of loose sand underlain by 
dense sand. Electrical strain gauges v;ere used to measure the deflec- 
tions along the length of the pile. Tests v;ere also conducted on 
partially embedded piles with an overhang of 8 feet, which showed 
considerable increase in lateral deflections at the top. 

Peck and Davisson (l962) reported the tests conducted on steel 
pile of 54 feet in length and embedded in silty soil, with thepjile 
tip resting on hard rock. The analysis of the data based on linearly 
varying soil modulus showed good agreement between predicted and 


observed values 
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Mori (1964) conducted full scale lateral load tests on steel 
pipe piles embedded in loose to dense sand. The deflections at the 
top v/ere measured under the action of lateral loads. The analysis based 
on the' assumption that the soil near the ground surface reached ultimate 
resistance showed that, the measured ralues were in aggrement with the 
expected values, 

Kerisel (I965) described the results of tests conducted on 
steel caisson embedded in clay. Values of displacements, slopes, 
moments and pressures were measured along the embedded length of the 
pile by using electrical strain gauges. The point of maximum bending 
moment was observed to be not at a greater depth from the top. Eeese 
(1965) analysing the above results showed that the analysis can be 
carried out on the assumption of linearly varying sail modulus. He 
also showed that, the approximate exp)ression given by him (Reese, 1958 ) 
to determine the ultimate soil resistance at the top is reasonably 
valid. 

Eubo (1965) reported the exhaustive results of the laboratory 
studios made on model steel piles. On the basis of the experimental 
results a non- linear soil pressure-deflection relationship which varies 
linearly with depth was proposed. It was also observed from the test 
results that the soil at the top goes into plastic region. 

Milaan and Eousberg (I965) suggested an optical method which 
can be adopted for measuring the horizontal pressures coming on piles 
when subjected to lateral loads. 
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Prakash and Saran (1967) performed model tests on aluminium 
piles in cohesive soils. The tests were conducted on single and on 
4 and 9“PiTe groups. It was concluded from the tests that the 
interference increg,sed as the spacing decreased in the direction of 
horizontal loads. 

Affad and Petrosovits (1969) presented the results of 
laboratory tests on model rigid piles embedded in cohesionless soils 
T^ith different batters. The test results show that the variation of 
coeffioient of soil modulus is considerably af footed by batter of piles,. 

Davisson and Salley (1970) presented a relationship between 
the behaviour of a single and group of laterally loaded piles from the 
tests conducted at Arkansas River Navigation Project, They observed 
for a single pile (a) a triangular variation of subgrade reaction for 
initial loading (b) at low loads a sharp decrease in subgrade modulus 
with increasing deflection, and not much change in soil modulus at 
higher loads, (c) 100 cycles of load caused on increase in deflection 
to 1.7 times that observed for first cycle, 

2.3 MALYTICiL INVESTIGATIONS 

In this section a comprehensive review of the analytical 
studies is presented. Most of the analytical studies reported in the 
literature treat the problem of flcsural behaviour of pile as a beam 
or beam-column on elastic foundation. Hetenyi (1946) and Vlasav and 
Leonte’v (1966) have presented the solutions for the fundamental 
differential equation for several boundary conditions and for variable 
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moment of inertia and soil modulus. 

Earlier works of Titze (l932)t Ratzersforfer (195^, vide 
Davisson, I96O) , Chang (l957) Pender (1947) s-H deal with the 
analysis of the flosural hehaviour problems on the basis of beam on 
elastic foundation concept , 

Palmer and Thompson (1948) have proposed the finite differenoe 
approach to the problem of laterally loaded pile based on soil modulus 
concept and have shown that this approach can be used with ease even 
for the most complicated variation of soil modulus. 

Subsequently, Barber ( 1955) » Palmer and Brown (1954)» Mason 
and Bishop (l954) and McClelland and Pocht (l958) presented numerical 
solutions by making use of finite-difference approach^ considering both 
constant and linearly varying soil moduli. 

In 1955» Terzaghi in his classical paper on subgrade modulus 
proposed that for preloaded clays, the soil modulus can be taken as 
invarient with depth whereas for sands it can be assumed as linearly 
varying with depth, 

Reese and Matlock (195^) proposed an iterative method of taking 
the soil modulus as the secant modulus in the non-linear range of the 
load deflection curve and make use of concept of beam on elastic 
foundation. For the solution they used the finite difference technique 
and presented the non»dimensional charts for deflection, moment and 
shear force values, for linearly varying soil modulus. 
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Davisson (i960) presented the solution for laterally and 
axially loaded pile TiTith linear varying soil modulus and with invariant 
moment of inertia and axial force with depth, hy electrical analog 
method, 

Matlock and Ingram (1965) have presented a general method, 
for solving the prohlem of "bean and beam-column on elastic foundation, 
based on finite-element mechanical models, 

Davisson and Gill (1965) presented electrical analog solution 
for laterally loeided pile in a two- layered soil system. The stiffness 
of iiie surface layer is defined in terms of the underlying layer. The 
study clearly broughtout quantitatively the overwhelming influence of 
upper few feet of soil on the flexural behaviour of laterally loaded 
piles, 

Mori (1964) presented solution for laterally loaded pile in 
elasto-plastic soil. Plastic resistance in the plastic region was 
assumed to be constant with depth. In the elastic zone th soil 
modulus was assumed to be constant with depth. 

Broms (1964 a, 1964 b) presented, the "ultimate resistance" 
analysis of laterally loaded piles in both oohosive and cohesionless 
soils. It is assumed for laterally loaded pile the plastic resistance 
develops all along the length of the pile for rigid pile. Ihexeas for 
the case of flexible pile, the ultimate load carrying capacity is 
arrived at from the consideration of failure of the pile because of 
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development of plastic resistance in the top region of the soil. 

Based on this previously published data were analysed and good 
agreement was seen iDetween the predicted and measured values. On the 
basis of comparison lath the field and laboratory measurements he 
also showed that the concept of soil modulus is reasonably valid at 
working loads, 

Spillars and Stoll (1964) presented analysis of laterally 
loaded pile by making use of Mindlin' s theory for a point load in semi- 
infinite continuum. In the analysis the effect of plastic region was 
also taken into consideration. 

• Poulos (1971 a) presented the solution for horizontal 
displacement and rotation of a vertical pile subjected to lateral 
loading and moment, and situated in a homogonious, isotropic and semi- 
infinite elastic mass using Mindlin' s theory approach, Numerical 
results are presented for a wide range of pile flexibilities and 
length to - diameter ratios, for both free -free and fixed-free casee. 
Comparisons between the elastic solutions and the corresponding 
solutions obtained from Winkler' s theory show that, the latter 
considerably overestimates the displacement and rotation of the pile, 
"tmt gives a very reasonable estimate of the moments in the pile. The 
elastic analysis is also extended to include the effect of local yield 
between the soil and the pile, Poulos (1971 b) also extended the 
analysis to group of piles. Based on the same approach recently, Poulos 
said Madhav (l97l) presented the results fcxrr axially loaded piles and 
laterally loaded piles. 



16 


Mailock and G-ru'b'bes ( 1965 ) analysed the problem of laterally 
loaded pile in elasto-plastic soil. The non-linear nature of the 
deflection curve is considered in the analysis. Extensive numerical 
results are presented to bring about the effect of various parameters, 

Reddy and Valsangkar (1966, 1967» 1968, 1969# 1970# 197'V) made 
a comprehensive study on the flexural behaviour laterally and laterally 
and axially loaded pile problem. For many oases solutions have been given 
in closed form, series solution and also by energy method. The effect 
of variation of soil modulus, variation of moment of ineirtia of the pii©# 
and variation of axial force has beon presented in these v/orks. 

Solutions are also presented for the problem of laterally loaded pile 
in both two-layored soil system and Elasto-plastic soil. 

Soy ( 1970 ) considered the different stages, starting from 
fully elastic to fully plastic soil behaviour, and presented the 
solutions for laterally loaded pile problem. 

Domenico Lalli (l97l) presented the solutions for the problem 
of pile with generalized loading in non-linear soils by large deflection 
matrix methods’ . The variation of cross-sectional properties of pile 
is considered and numerical results are presented using the proposed 
method for a particular case. 

Basudhar (1971) obtained solutions for the laterally and 
axially loaded pile problem in series form and the results are 
presented in non-dimensional forms. The effect of variation of 
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axial force and soil modulus is considered and exiiaustiye nujnerical 
results are presented for layered soil, elasto-plastic soil and for a 
pile driven throu^ a compressible strata. 

2.4 GEIIEEAL CONCLUSIONS 

Based on the revieu- of the experimental and analytical 
investigations certain generalized ponolusions are drarJn he low, 

(i) The flexural behaviour of a pile is gewomed hy the 
properties of .top region of the SDil of about 10 to ,15 diameters of 
the pile, 

(ii) Coefficient of soil modulus and plastic resistance are 
considerably influenced by different factors, viz ; Stiffness of the soil, 
remoulding, repetition of lateral load, vibrations du3fLng giving, 

consc lidation coefficient and size of loaded area. 

(iii) For a more realistic analysis of pile problem, it is 
necessary to consider both elastic and plastic state of the soil. 
Especially at hi^er loads it is known that the soil at the top goes 
into plastic region. 

(iv) When in addition to lateral loads, axial load also is 
acting, the flexural capacity of the pilo considerably decreases. 

(v) The soil modulus concept is reasonably valid at working 

loads.' 

(vi) Boundary conditions play an important role in the 


flexural behaviour of pile 
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(vii) Field tests have ' shown that Terzagni' s (1955) 
reoomniendations are fairly valid for soil modulus variations, 

(viii) Brom' s approach oan he used for finding the 
ultimate carrying capacity of the pile with reasonable accuracy. . 

(ix) As the review shows that the behaviour is considerably 
influenced by a variety of parameters and as such a probabilistic 
rather than deterministic formulation will be an ideal one. 

(x) Prom the tests of Poulos (1971 a) it can be concluded 

I 

that the winker’s approach can be tised with reasonably good accuracy 
to obtain the bending moment with depth for a pile problam, and the 
elasticity theory has to be used to get a fairly good estimate of 
deflection and slope. 





cHipm 5 


GEHERILIZSD SOLUTION FOR LATERALLY MD AXIALLY LOADED TAPERED- PILE 

IIT LAYERED SOIL 


5.1 GEEERAL 

The necessity of flexural analysis of axially and laterally 
loaded piles arises from the use of these piles in marine structures, 
oil drilling rigs, cranes and other similar structures wherein these 
piles have to resist large lateral and axial forces. It is well known 
that, when in addition to axial .forces, the pile foundations have to 
resist large lateral forces, the flexural capacity of the pile foundar- 
tions oonsiderahly decreases, because of "Beam*coluinn effeat ” , Field 
investigations of Arup (1944), Evans (1955)» Barber (>955), Mcnulty 
(’'956)) Peck and Ireland (1961) and Francis et.al. (1961) support 1his 
view, Hetenyi (1946) bas presented many analytical solutions for a 
beam-column resting on elastic foundation for different boundary 
conditions, 

Eventhcfugh the overwhelming effect of the upper layer on the 
flexural behaviour of piles was known longback, analytical studies were 
presented only recently by Davisson and Gill (1963) for the case of 
uniform cross-section pile in a layered soil Trd.th constant soil 
modulus in each layer. Recently Reddy and Valsangkar ( 1968) considered 



the general variation of soil modulus in the fom of a second degree 
polynomial and presented the solution for a pile of uniform cross-- 
section emhedded in a layered soil and sahjeoted to lateral load only. 
Basudhar (1971) obtained the generalized solution for a pile of uniform 
cross-sectioUj embedded in a layered soil, acted upon by a lateral load, 
moment and an axial load, considering the effect of variation of axial 
force along the embedded length of the pile. 

In many instances step-tapered, uniformly tapered and stepped 
piles are used for transfer of load from superstructure to subsoil. In 
case of step**tapered piles, because of corrugated surface and the various 
steps, the load carried by the skin-friction component is considerable 
and as such these piles are used as friction piles (D’ Appolonia and 
Hribar, 1963) . Experimental studies by Mcnulty ( 1956) and. Peck and 
Ireland (1961) on Raymond standard and step-tapered piles show that 
reinfo3X!Gd piles have better resistance to lateral loads when compared 
to unreinforced piles. 

Review of the literature also reveals that, the analytical 
studies taking into account variation of moment of inertia are very few. 
Reese and Gingbarg (1958) lave proposed finite difference method for the 
analysis of laterally loaded piles with abrupt changes in moment of 
inertia and linearly varying soil modulus. 

As the available literature confines to only uniform cross- 
sectioned piles in layered soil, analysis has been made in this chapter 
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to Ecrrive at a generalized solution for the case of a tapered pile 
emhedded in a layered soil and acted upon hy generalized loading. 

5,2 STATEIOTT OF IBE PEOBiai 

Fig. 3.1 (a) shows a uniformly tapered pile with circular 
cross-section whose moment of inertia varies with fourth poT/er of x. 

Fig. 5.1 (t) shows a free head tapered pile T/hicli is acted upon hy a 
lateral load, moment, M^, and an axial load, P, at the top. 

Fig. 5*1 (c) shows schematically a fixed head pile, wherein it is assumed 
that the pile is free to translate and is restrained against rotation. 

Fig, 5*2 (a) and'Fig. 3,2 (h) show the general representation of 
laterally loaded pile problem along with the possible complexities in 
the behaviour of the soil. In Fig. 5*2 (a) the form of springs, indicate, 
that the stiffness of the soil varies with depth and deflection of the 
pile. The St. Venant' s body elements at the top indicate that there is 
some limit to the soil reaction which can be developed in the top region 
and the Newtonian dash-pots represent the time dependent viscous effects 
of the soil. It should also be noted that the repetition of load alters 
the behaviour of the soil pile system. However taking all these factors 
into account makes the problem involved mathematically, as such the 
analysis presented herein takes into account elastic nature of the soil, 
variation of moment of inertia and variation of axial force along the 
embedded length of the pile. Fig. 3.2 (b) shows the probable variation 
of soil reaction along the embedded length of the pile. 
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3*5 (a-) and Pig. 3*5 (^) depict the linear and non-linear 
soil pressure-*def lection relationships respectively. Even if the p-y 

relationship is linear may vary t/ith depth and , hence a general 

s 

relationship would he, p = E (x) y, If p-y relationship is not 

s 

linear as shoT/n in Pig. 5*3 (h) , thcgi the soil modulus will he a function 
of hoth X and y. Matlock and Roose (1956) suggested the use of the 
secant modulus and repeated application of the clastic theo3^y to arrive 
at the solution of this non-linear problem. 

Since the concept of soil modulus is an important factor in the 
following chapters a brief discussion of the same is given below. 

Generally, for all soils, the value of coefficient of soil 

modulus will be smaller at ground surface than at some depth below the 

surface. According to Terzaghi (l955) the soil modulus is constant with 

depth for pre-loaded clays as shown in Pig. 3*4 (a). In many cases there 

may exist at the top a stiff crust due to dessication or some other 

mechanism. Then the probable variation of soil modulus is as shown in 

Pig. 3*4 (t) . Pig. 3*4 (c) represents the probable variation of soil 

modulus for the case of a weak soil layer over a hard layer (Davisson 

and Gill 1963)* Po^? sands E is assumed to vary linearly with depth as 

s 

shown in Pig. 5 *4 (d) . Por normally loaded clays also variation of E 

s 

is given by Pig, 3,4 (d) (Davisson, 196 O) . 

Pig* 3*5 Pig, 3*6 show a free-head pile and a fixed -head 
pile respectively, in a two layered soil with second older polynomial 
variation of soil modulus in each layer, which can be considered to 
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be the most general variation (Reese and Ilatlock, 1956) . However, in 
the present work it is assamed that the soil modulus is constant in 
each layer, so as to have mathematically elegent solutions. As stated 
earlier, this case corresponds to the variations as depicted in Pigs, 

3*4 (a), (b) and (c) . 

As in the present case, the projected width of the pile 
varies with x, a brief discussion regarding the effect of width on 
coeffioiont of soil modulus is given. 

Torzaghi (l 955 ) has presented a detailed disciiasion regarding 

the effect of width on the coefficient of soil modulus. Pig. 3*7 shows 

a pile embedded in cohesive and cohesionless soils alongwith the 

presaire bulbs for piles of widths B and nB. Por linear behaviour of 

soil and the same intensity of loading, the pressure bulb for a pile 

is 

of width B for both clays and sands/as shown in Pig. 5«7 • Because of 
this the deflections for a pile of width nB will be n times the 
deflections for pile of width B. Thus the coefficient of soil modulus 
for a pile of width B is equal to the coefficient of soil modulus for a 
unit area multiplied by width B. Because of this fact, the soil pressure, 
p, par unit length of the pile remains invariant with depth (Davisson, 
i 960 ). However, the above reasoning is only valid for elastic behaviour 
of the soil and the effect of width is different in the non-linear 
range of p-y relationship. 

The laboratory studies of Matlock and Ripperger (1957 B) 
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confirm Terzaglii’ s recommendations regarding the effect of lyidth on 
the coefficient soil modulus. They also concluded that the deflections 
of a pile for a given lateral load should he the same regardless of 
pile diameter. Hence in the present chapter, taking the validity of 
Terzaghi* s (l955) recommendations, ana.lytical solutions are obtained 
for constant soil modulus. It should be noted, that the analysis 
presented herein should be used only for linear p-y relationship. 

The variation of moment of inertia with depth x for the 
tapered pile shown in Pig. 5*8 (a) is given by ; 

l(x) = Iq(i + bx)^ ....(5.2.1) 

0*" i 

in which = moment of inertia at the top 5 b 
L r: length of the pile, 

and 0 = ratio of bottom diameter to top diameter. 

Even though in the present analysis a continuous variation of moment of 
inertia is assumed, the same analysis can be used for finding the 
flexural response of stepped and step-tapered piles, as in general the 
steps at each section are small and one can appro 2 d.mate the pile by one 
which is uniformly tapered. 

The variation of axial load depends on the distribution of 
skin ftiction which in- turn depends on the properties of the soil, 
pile etc. For clays it is a time dependent phenomenon whereas for 
sands it is invariant with time. 
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It should he noted that similar to contact pressure distribution 
at the interface the aodLal force distribution is governed by many factors 
and as such generalization regarding the distribution is not possible. 

As such in this trork it is assumed to vary parabolically as shown in 
Fig. 5*8 (t) . Thus the axial force at any depth x is given by ; 

P(x) = P(l+bx)^ .... (5.2.2) 

in which 


P a axial force at the top of the pile. 

Therefore the axial force at the tip is given by ; 

P(x)'= Pc^ .... (5.2.5) 

at 35 sL 

So the load transferred due to skin friction for this type of variation 
of axial force is equal to P(l-o ). For c = 0.6 (usual value) the load 
transferred due to sldn friction is 64 percent of total load coming at 
top. This particular form has been chosen as this would lead to a 
mathematically simple analysis . However, it should be noted that, at 
axial loads corresponding to O.5 to O.5 the effect of variation of 
axial force on the fl'>xural behaviour is practically negligible (Reddy 
and Valsangksr,1968) 0.nd. as such in these range of axial loads it is 
immaterial which form of axial force is assumed. Thus the assumed 
variation thou^ being particular, the results can be used for the load 
ranges upto O.5 without much errors. 

5.5 MALYSIS 

The governing differential equations for the uniformly tapered 
pile shown in Fig, 5«8(a) wherein, constant soil modulus is assumed in 
each layer are given by (Hetenyi, 1946) : 
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dx 


• dx 


El(x) 


El(x) 


d^y 

dx 

+ dT 

dy 

L J 

d^ 7 { 
dx^ J 

1 

d 

dx 

' 3 -yn” 

— -J 


+ K = 0 


+ = 0 


O^x^^L. (3.5.1) 


L^:^x4 L 


. (3.5.2) 


It is olrrious from equations ( 3 * 3 . 1 ) and (3.3.2) that, even 
for' the simplest ease- when' 'the axia,l force remains constant with depth, 
it is difficult to obtain the analytical solution. However, for the 
particular axial force variation shown in Pig. 3.8, it will he ^own 
herein that, it is possible to arrive at a closed fonmi solution. 

Substituting for l(x) and P(x) from equations (3.2.I) and 
(5.2.2) into equations (5.3.l) and (3,3,2), the above equations 
transform to ; 


dx 


[EI^(l+bx) 


;,2 

d y 

I?’ 


d 

dx 


[P( 1+bx)' 


dx 


3 + ^^o^u 


0 


(3.3.5) 


01 X <L^ 


2 ^ 2 

Afelo ^ ^ p( 1+l^2c) 

dx , 2 Jj dx I— 

dx 


)2 dy ^ 

dri+Vi = ° 


.... (3.3.4) 


^ii ^ 1^ 


To feoilitate the ntunerical computations, equations (3. 3 . 5 ) ( 5 . 5 * 4 ) 

are trensfermod into non-dimensional form by making use of the 
following non-dimensional parameters ; 
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Depth coefficient Z = ^ 

Jl 



and masiimim depth coefficient! Z* = 

Jl 

in which R = relative stiffness factor 


• » • • (3*3 * 5 ) 
.... ( 3 . 5 . 6 ) 

• (5.3.7) 


In this work, it is taken as ; 

R = *••• (5.3.8) 

In. the present analysis, since elastic theory is used, "by 
assuming the validity of the principle of superposition the actions 
of (i) axial load, P and lateral load, (ii) axial load, P and 
moment, are considered separately for free head pile. Thus if y^ 

and yg are deflections due to (i) and (ii) respectively, then the total 
deflection is given hy : 


y =« + ^2 


.... (5.3.9) 


Similarly, the bending moment, M under the combined action of (i) and 
(ii) is ^ven by j 


M = + Mg 

in which = B.M. due to (i) 

Mg = B.H. due to (ii) . 


(5.3.10) 
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dimensional deflection and moment coefficients axe defined as follows. 
Case (i) Fxee-fxee end conditions pile, only Q and P acting at the 

"U 

top ; 


t ^1®^o 

Deflection coefficient, Y. = ^ 


.... (5.5.11) 


and Moment oo efficient, 


M 

t 


• • • • (5.5.12) 


in which and axe dimensional deflection and tending momcait values 
xespectively. 

Case (ii) Free -free end conditions pile, only 11^ and P acting at 
the top ; 


, y2Eio 

Deflection coefficient, Y 2 = g v .... (5.5*15) 


V 


and 


, “2 

Moment coefficient, ]yi2= 

t 


.... (5.5.14) 


Case (iii), Fixed-fxee end conditions pile, with and P acting 
at the top i 


i y-ifEio 

Deflection coefficient, Y.^ = ^ ; 


QtE^ 


“if^^o 

and Moment coefficient, M. „ = Tyrr — 

ii 


.... (5.5.15) 


.... (5.5.16) 


y^^ = dimensional deflection and M^^ = dimensional bending 


in which 
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moment for fixed - free pile. 

By making use of the above non-dimensional coefficients 
equations ( 5 * 3 * 5 ) siad ( 5 « 3 » 4 ) reduce for the above case (i) to : 


^2- [ (1+d’Z)^ 3 ^ ^ f (l+d'Z)^ = 0 

dZ dZ 


dZ " ‘ ^o^ 1 u 

oiz<z^ .... (3.3.17) 


and 

,2 


2 t 

.2 . d Y . dYi 

-2 t(n-a.z)'^^ ]+!!-[ (i+d'z) ....(3.5.18) 


1 ^' 


in uhioh 


dt = 


0-1 

(WT 


u 


m 

E^o 

Kj 


^..( 5 . 5 . 19 ) 


and 


r 

o El 

0 


Since the above differential equations are of Euler-Couchy type these 
are transformed into equations with constant coefficients by the 
following sabstitutions, (ince, 1956, Kamke 1959, Kosko 1965) *. 

(l+d'z) .... (5.5.20) 

aM \ ....(3.5.21) 

Then equations ( 3 . 5 . 17 ) and (5.3. is) transform to ; 
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^ + 2 + ( 4 . - 1 ) — ^ +(,(.'-2) 


dY‘ 


dv 


3 


dv 


In • 

+ r AY, = 0 
dv o 1 u 


.... (5.3.22) 


and 

dV dV dV , 

— ^ + 2 —^ +(ijj -1)--^ +(il; -2)— I- 4- ^Y =0 .... ( 3 . 3 . 23 ) 

dv dv''^ ’ dv dv 


t U 

in whiith 4 > = — f— ^ 

(a )^ 


and 


. » • • ( 3 . 3 . 24 ) 


(d' )' 


4 


Eofuations (3.3*22) and (3.3.23) a^e homogeneous ordinary linear 
differential equations irith constant coefficients for which closed 
form solutions can he arrived at as follows ; 

Writing in the operator form we have 


. [A2d\ (if»'-l)D^+(i|^'-2)D + r 


■Xl'Y, ='0 
1 u 


.... ( 5 . 3 . 25 ) 


and[D'^+ 2 l \ ( ,f,*-l)l)\( V -2)D + X]y1j3_ = 0 


.... (3.3.26) 


The solution of equations ( 3 . 3 . 25 ) and ( 3 . 5 . 26 ) are assumed to “be 
of the form s 


I 



a V 

CG 


.... ( 3 . 5 . 27 ) 


Then the characteristic equations are obtained as • 
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CJ + 2(3 + ('>|^'“l)c[+(v "*2) oi H- A=0 ...,( 5*3 *28 ) 

and 

/+ 2a^ + (ii''-l)cx^+ (V-2)a + X = 0 ....(3.3.29) 

Equations (3.3.28) and (3.3.29) are sol-ved by Brown's metbod for 
solving Eourth order polynomial (Mc'Gormio and Salvodari, I964) and the 
roots have been chocked by B air stow ’ s method (Csrnahan, Luther and. 
Wilkes, 1969) . 

The general form the roots to is ; 

Oj to =» — •^+0 + iT- .... (3.3.3‘^) 

in which numerical values of cr and x depend on the values of 1^* 
and X . However in all cases the form of roots remains the same as 
equation (3.3.30). 

The general solutions of equations (5.5.22) and (3..5«25) after 
re substituting the values of v and n are (ince 1956, Kamko 1959> 
Kosko 1965) : 

t 

Y^u= (l+^'2) A Cosh[aln(l+d'Z) ] Cos [t ln(l+d«Z)] 

+B Cosh[u ln(l+d'Z) ] Sin[ Tln(l+d'z)] + C Sinh [uln(l+d*z)] 
Cos[T ]n(l+d'Z)]+ D Sinh [aln(l+d»Z)] Sin -fTixL( i+d» Z)] } 

.... (3.5.31) 
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and 

f ! ! I 

= (l+a'Z) 2{A Cosli[o>ln(l+d Z)]Cos[t' ln( 1+d Z) ] 

+ B Cosh[CT ln(l+d Z)]Sin[T In (l+d z)] + C Sinh.[ CT*ln(l+d Z)] 


Cos[t ln( 1+d Z)] +D Sinh. [ <^ln(l+d Z) ] Sin[T ln( 1+d Z) ] > 

• • • • (3 •5»32) 

1) t t I 

The constants A, B, C, D, A , B, C and D are evaluated from the 
■boundary conditions at the tip and top of the pile and from the ' 
compatibility conditions at interface between the tr/o soil layers. 

For the case (i) i.e, for freehead pile acted upon by daear force, Q,^, 
and axial load, P, the boundary conditions are j 

I. at xs=0 i.e, at Z=0 

.... (5 •3.*553-} 

(3.3.53 ) 

U J kJ oJ.’ ~ 

<3-Vu. 

^ [-El(x) — ^ ^ .... ( 3 . 5.34 a) 

dx 


a) Bending moment = 0 

El(x) » 0 

dx^ • 


or 

and 

- \ « -r-, 


2 * 

ill- 

dZ^ 


=• 0 


3 t 2 « 

^ ^1U 

+ 43 + u 


dZ 


dZ 


dY 


J 

1u 


or 


dZ 


1 


....( 5 . 5 . 54 ) 
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II at X = or Z =s . The 4 oompatihility conditions 
a) Hof lection oompatihility. 


I » 



h) Slope oompatihility 


dT 


In 


dY 


11 


dZ 


dZ 


c) Bending moment oompatihility. 


2 t 

d Y 


1n 


Vi 

dY 


11 ' 


dZ 


dZ 


d) Shear force oompatihility. 


and 


3 ’ 3 ' 

d-^Y. d^Y_ 

In 11 


dZ 


3 . 


dZ 


3 " 


III, at X = L ox Z = Z 


a) B = 0 


2 * 
d Y 


X .e. 


11 


= 0 


dZ 


and 

h) S.P = 0 


5 ’ 

d^Y^^ 4d 
+ — 


2 ' 
d Y 


11 


U 


+ 


dZ^ (l+d'z’) dZ^ (l+d'z')^ 


fill 

dZ 


* 0 


axe t 


.. (5.3.35) 


. (3.3-.36) 


.. (5.5.57) 


... (3.5.38) 


... (5.5.39) 


... (5.340) 
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By usiiig tlie aBove 8 TDOundary conditions, equations (3*3*33) 

I r I t 

to ■fctie constants A, B, G, B, A,S,C and I) are evaluated. 

For the case (ii) i.e, free head pile acted upon hy a moment, 
and an axial load, P at the top and case (iii) i.e, fixed head pile 
acted tip on hy a lateraA load, Q,,^, and axial load, P, the general form 
of the solutions, equations (3*3«31) and (3»3*32) remain same. The 
Boundary conditions at x = L and compadihility conditions at x = 
also remain same. However the Boundary conditions at x = 0 are ; 

Case (ii) For free head pile, and P acting at top : 


at X = 0 
a) S.P.= 0 


I.e, 


5 ' 2 ' ' 

d^Y„ , d Y„ dY„ 

2u . . T ’ 2 u , tt 2u 
— = — +4d — 5 — + H 


dZ' 


dZ 


dZ 


= 0 


• • •,(3»5 *4 ^ 

....( 3 . 3 . 41 ) 


B) B.1I = .... ( 3 . 3 . 42 a) 

2 » 

d Y 

i« 6 . o “ 1 .... (3»3»42) 

dZ 


Case (iii) for fixed head pile, only Q,^ and P acting at top : 
at X = 0 


a) slope a 0 


dY 


1 ,e. 


1f 


dZ 


= 0 


• • • • (5 •5»433') 

(3.545) 
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.... ( 5 * 3 » 44 a) 

.... (3.544) 

3.4. EESULTS MD DISCUSSIOH 


■b) S.P = Q, 


31 2 » ' 


Based on the generalized solutions presented above, 
e:5diaustive numerical results are obtained for the three different 
cases viz : 

1) Brec-Pree pile with only P md Q,^ acting at top 

2) Free -Free pile v;ith only P and acting at top 

and 3) Pixed-Pree pile with P and acting at top. 


Matlock and Reese (196O), Davisson (1960), Davisson and 
Gill (1963)^ Reddy and Talsangkar (1968) have shown that for all 

I 

practical purposes, the laterally loaded pile with Z = 2 or less ‘ 

I 

behaves as a rigid member and with Z = 4 or more behaves as a 

flexible member of infinite length fox both free-free and fixed-free 

end conditions. Hence for obtaining the numerical results in this 

t 

work, Z is taken as equal to which is a suitable approximation 

for most of the piles. Also, the value of 0 is taken as 0.6, which 

is the usual value for most of the field cases. For numerical 

computations the axial load ooefficient IT defined by equations 
(3*5»'19) is expressed as a function of the critical load coefficient 
of a free-free pile of uniform cross-section without skin friction* 
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Results are presented for the value of U = O.5 and 0,6 TJ^^,T7hore 

TJ = 0.9195f is the Buckling load coefficient fox a pile of 

length Z*= 4-0* in the analysis it io tacitly asooncd that small 

deflection theory for hoam-column is valid till this range. However, 

it should he noted that , since the solution obtained viz ; Equations 

(5 *3 *51) and (7.5*32) , are general in nature, any value of Z’,c and U 

can he used for computation purposes. Eshaustive results are obtained 

for r^ = 0,5, 1.0, 2.0, 4*0;aiid 6.0 which represent a wide range of r^ 

found in the field. For each value of r results are obtained for 

0 

different values of Z^ = 0.2, O.4 and 0.6. 

The first set of results were obtained for r^ = 1 corresponding 
to uniform soil and the results W'-'re found to be in ■conformity with 
the results of Valsangkar (1969)* 

Figs, 5«9 snd 3*10 show the offoct of axial force on deflection 
and moment coefficient values for free-free pile for 0.2 and r^=2,0. 
It is seen from the Pigs. 3.9 and 3. 10 that there is a considerable 
increase in deflection and moment coefficient values with increase in 
axial force. Also as the axial load approaches the buckling load, the 
deflection and moment coefficient values get considerably magnified. 

Pig, 3‘11,gi''^®s similar information for fixed free pile, Por tliis case 
also for numerical computation H is assumed to be equal to 0,9193« 

_ increases from 0.3 to 
0.6, the increases in deflection and moment coefficient values for 

free-free pile are ; 47*5 and 69,6 per cent for maximm 

1 

and respectively, and 65.7 percent and 29. 0 per cent for maximum 

Y2 * and M2' respectively. The corresponding percent increases for 

fixed-head pile are 7 per cent and 5 percent for maximum Y^^ and 
respectively. Thus it is found that, under 


Prom these results it is seen that, as TJ/U 
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the same axial loads the magnification in deflection and moment 
coefficient values in case of fixed-head pile is considerably lesser 
than tha.t of free-head pile. 

Figs. 3*‘12 to 3*17 represent the effect of r^ on deflection 

11(11 I 

and moment coefficients, and J.h , and i'L and M-jf and Y, „ 

. ' 1 "I ' 2 1f 

for = 0,2 and = — = 0.5 and 0.$ . It is scon in general from these 

^cr 

graphs that there is considerable increase in deflection and moment 
coefficient values for r^^ less than unity and there is considerable 
dampening of deflection and moment coefficient values for r^ values 
greater than unity. 

Figs, 3. 18 and 3,19 depict the effect of varying r^ on 

( t 

deflection coefficients, Y^ and Y2, at Z = 0 for Z^=0,2, O.4 and 0.6, 

IT 

for ~ = 0.3 and 0,6. Here again it is seen that for values of r 
or 

less than unity there is considerable increase in deflection coefficients 
with increasing values of and for values of r^ greater than unii:^ 
there is oonsidorablo decrease in deflection coefficients with increas- 
ing values of Z^ . llso it is seen that for values of r greater than 
° . 1 o . 

unity, there is relatively more docrease in deflection coefficients, 

I t 

and Y^vZ^ changing from 0.2 to 0,4 than when Z^ changes from O.4 

to 0,6. From this it can bo reasoned out that for r^ greater than 

! 

unity as Z^ approaches Z , the rate of decrease in deflection 

I ( 

coefficients, Y^ and Y^ considerably decreases, Ihis brings about the 
beneficial effect of having a hard layer at the top upto a certain 
depth, to control the deflection and moment coefficients. 



Pigs. 5*20 and 5 *21 show the effect of on majdimiin and 
for d/lJ K 0,5 and 0,6 , for different values of Z = 0,2y0,4 
and 0,6 . Hero again for r values less than unity the maxiimm 

) t 

moment coefficients, and increase v/ith increasing Z^ , For r^ 

t 

values greater than unity upto 4»0, the maedmum decreases with 
increasing Z^, whereas for r^ values greater than 4*0 there is no 

I t 

exact correlation hetiTeen r^ and and for .• .different values 

of Z.J. This may he reasoned out from the fact that the maxinum moment 
does not occur at a particular depth for all the values of Z^ . In 

t 

Pig. 3 *21 it is scon that, the effect of Max. 1/!^ is similar to 

! 

that on M, . 

1 

Pigs. 3 *22 and 5»23 show the effect of increasing r^ on 

I I 

deflection and moment coefficients Y.^^ and at Z = 0, for 

U 

=jr— = 0,5 and 0,6. Here" also as expected for r values 3e ss than unity 
^cr ° 

the deflection and moment coefficients increase with increasing Z^ and 
for 2 ^ vahioG greater than unity, decrease with increasing Z^. 

3.5 CONCniSIOMS. 

Based on the foregoing analysis and the numerical results 
presented, certain genora^lized conclusions can he drawn. The problem 
of tapered pile under generalized loading in layered soil presents the 
most general case of the pile problem met with; since the case has 
varying pile section properties, soil properties, and axial force. 

The analytical solutions obtained being in closed form can he used 
as a check on the convergence and numerical accuracies, obtained 
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CHAPTER 4 


GENERALIZED SOLUTION EOR AXIALLY AMD LATERALLY LOADED TAPEERED PILES 

IN ELASTO-PLISTIG SOIL 

4 . 1 GENERAL 

In Chapter 5 the studies were hased on elastic theory which 
is valid at working lo^s, wherein p-y relationship was assumed to he 
linear according to Broms (I 964 a »1964 h) . ■ However, at higher loads 
the soil hohavGs non-linearly as reported hy ; ^agner 1955> ^'itzgihhon 
1956» Matlock and Eipperger 1957 an<t Mori 1964* To account for this 
in the present chapter the soil is assumed to he represented hy linear 
comhination of St. Venant' s body and linear spring as shown in Eig,5.2. 
Pig. 4*2 represents the p-y relationship for this model. Mori ( 1964 ) 
has presented the solutions for a laterally loaded pile problem with the 
above assumption for infinitely long pile and the correlation with field 
data was scorn to he reasonably good. Earlier Starzewski (1959) observed 
from the experiments conducted on beams loaded at edge and resting on 
soil that, the measured soil pressures vjorc in bet./oon the pressure 
predicted by soil modulus concept and that obtained from assuming the 
p-y relationship as shown in Pig. 4*2. 

Evonthough the fact that the soil at the top region of the 
pile goes into plastic region was recognized as early as 1955 
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(Wagner, 1953) > only recently Broms (1964a, 19646) presented the 

solutions for tho ultimate carrying capacity of laterally loaded piles 

taking into account the plastic resistance of tho soil. In tho analysis 

> 

it was assumed that either pile fails "because of excessive Bending 
moments or the soil goes into plastic resistance along the em'bedded 
length of the pile, Ho'wever in most of the cases usually only a fevf 
feet of the soil in the top goes into plastic resistance, which may not 
"be sufficient for the pile failure to take place, 

Reddy and Talsangkar (1970) presented the generalized solutions 
for laterally loaded pile in elasto-plastic soil. M.B, Roy (1970) 
sented studies on laterally loaded pile in elasto-plastic soil for 
constant soil modulus, and constant plastic resistance. Basudhar ( 1971 ) 
studied the effect of axial load also in addition to lateral load on 
the pile in elasto-plastic soil. Recently Domenico lalli ( 1971 ) 
presented the solution for axially and laterally loaded pile of varia’ble 
cross-section "by large deflection matrix method v/herein "bi-linear p-y 
relationship is assumed. The results are presented in dimensional form 
for a pile of particular cross-section. 

As, except for Lalli 's (l97l) 'i’'fork all other studies confine to 
uniform cross-section piles. In this study generalized solutions are 
presented for axially and laterally loaded tapered pile in elasto-plastic 
soil. The assumed variation of axial force along the length of the pile 
is shoTJn in Rig. 4»l(6). This particular variation has "been taken in the 
the analj^sis for tho reasons already discussed in Chapter 5« 
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4.2 STATEMENT OF KIE PROBLEM 

m 

Fig. 4*1 shows a circular uniformly tapered pile whose moment 
of inertia variation along the embedded length of the pile is ; 


I(x) = I^d+hx)"’- 


.... ( 4 . 2 . 1 ) 


in which 


^0 = 


nx). 


h = 


°rJ: 

' L 


c « ratio of the bottom diameter to the top diameter 5 
and = diameter of the pile at the top. 

In this schematic diagram the pile is acted upon at the top 
by (i) a lateral load and an axial load P. 

(ii) a moment and an axial load P. 


The above two cases have been considered separately in the present 
chapter. 

The variation of acdal force shown in Fig, 4*1 (l>) is given by; 


2 

P(x) a P(l+bx) 


(4.2.2) 


where P = axial force at the top of the inle. 

Because of relatively wealc nature of the soil and large 
deflections at the top, upto a certain depth the soil goes into plastic 
region as is concluded at the end of Chapbor 2. To account for this 
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plastic resistance offered by the soil the load deflection relation- 
ship is assumed as shown in Pig.. 4.2 in this analysis. Broms (1964 b) 
presented an analysis for the ultimate load carrying capacity of a 
laterally loaded pile in cohonionloss soil rath the variation of 
plastic resistance in the plastic zone with depth as 


^ult 


3B-YxKp 


.. .. (4.2.3) 


Reddy and Yalsangkar (l970) presented a solution for ultimate soil 
resistance for a laterally loaded pile embedded in cohesionless soil 
with the assumption that, a triangular wedge moves upward at the time 
of failure as in the case of cohesive soil, the corresponding variation 
of ultimate resistance is parabolic in nature, Basudhar (l97l) on the 
basis of the same assumptions as above presented generalized solution 
for axially and laterally loaded piles in .elasto-i^lastic soil. 


4.3 MALYSIS,, 

The governing differential oq.uations for the uniformly 
tapered pile shown in Pig. 4 .1? wherein it is assumed that, the bottom 
layer has constant soil modulus and the top layer goes into plastic 
region, are given by : 


2 d^ 


dx 


dx 


dx 


» 0 • 


• ••• (4*3*'l) 


- ^ 

d y^ 




dx 


[ElW ^ 53^1+ V- 0. I-xlL 


dx 


(4.5.2) 
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in whioh t = cocffioient of ultimate soil resistance. 


Substituting for l(x) and P(x) from equations 4 * 2.1 and 4 * 2.2 
into equations 4*3*1 and 4*3*2, leads to ; 


dx" 


^ [El^o+bx)"^ ^+ [P(i+bx)^ ]+ t = 0 .... ( 4 * 3 * 3 ) 

2 dx 


dx 


0 <x <L, 


1 


dx 


[EI^ 


ijl+bx)"^ -4V ^ [P(l+bx)^ ^1 ]+ = 0 .... ( 4.3.4) 

dx dx 


L.< X < L 
1 - — 


Igain to facilitate the numerical computations equations 4*3*3 
and 4*3*4 a^re transformed to non-dimensional form by making use of tlxe 
following non-dimensional parameters s 


Depth coefficient? Z = yi/R .... ( 4 * 5 * 5 ) 

Z^= •*.. ( 4 * 3*0 

and Minimum depth coefficient Z = L/R (4.5.7) 


in which R *= relative stiffness factor given by 




( 4 * 3 * 8 ) 


For' the oases of free-free pile and fixed-free pile, the non- 
dimensional deflection and moment coefficients are defined as follows; 
, ■ MqEIo 


yi = 


**** (4.3.9) 


tR^ 


M 


• • • • 


(4.5.10) 


'!-5 


t 



tR 


in v^hioh 

» 

« dimensional deflection 5 

and M = dimensional bending moment. 

By making use of the above non-dimensional coefficients 
equations 4.3.3 and 4.3.4 reduce to s 


2 

,2 /I i yl j o dy' 

i- [(1+d'z)'^ ti [(i 4 d'zf ^^> 1 ] + 1 . 0 

dz'^ dZ^ ■ 

0 <z< z^ 


.... (4.5.11) 


4 [(idd.z4^’-Hnfj[(i.d.z)2^M. y;,.o 


* 

d V 




Z^^ Z 


.... (4 »5 • 12) 


in whioli 


d> = 


B = 


0-1 

(L^) 

PE^ 

EI„ 


.... (4 * 3 . 1 5) 


The above equations 4 . 3*11 snd 4 * 3*12 arc of Eulec Cauchy 
typo and honoe these are transformed into equations with -“eonst ant 
coefficients by tho follovdng substitutions (inco 1956> Kamke 1959 
and Kosko 1965) ' 

•n 


1+d'Z = 


• » « • 


(4.3.14) 
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and 


V 


® ..t. (4.5.15) 

The final form of equation 4.5.11 is obtained as s 


4 

cL yi 


z 2 

d'^y* d y* dy* 

+ 2-yJ-^ +(i|;t-l) — +(^.-2)^ + X=0 
dv^ dT^ dv 


1u 


(4.3.16) 


in whioh 


TJ 


“ (d’.)^ 


^ (dt)^ 


• ••' ( 4 •3*17) 


X V 

SalDstitaixng « y !,^ + 4 - *5 *16 leads to ; 


4 3 2 ^ 

dT!, crY» d dYL^ 

-y " = 7T - if" 

dy dy dv 


.... (4.3.18) 


Equation 4.5*18 is a homogeneous ordinary linear differential 
equation v/ith constant coefficients for ivhich closed form solution can 
he orriyed at, as follov/s ; 

Substituting D = ~ 

[D^ + 2D^-h Gji'-l) dV (i/j'-E) E = 0 

The solution is assumed to be of the form* 


.... ( 4 . 3 . 19 ) 
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.... (4.3.20) 


Then .the characteristic eq,uation is given hy 
Y^''+ 2 Y^+ ('^'-1)y^ + Y= 0 


.. .. ( 4 . 5 . 21 ; 


The 4 roots of the equation 4.3.21 are y^ "fco 


’' 1-0 


Y 2«-1 


'-1~ y^~4i] 


.... ( 4 . 3 . 22 ) 


and the general solution is given by s 

ov -V (4 3 23) 

Yt = Afi 4. -Re J. no 2 + De^ 2 ^ 


yt^ = Ae + Be + Ce 2 


Putting a = ‘-h and 3 . = 


^ 4’J ' ' -9 
2 


Y Y — — ~3. 
y '4 “ a i 2 


a ± 1 


.... ( 4 . 3 . 24 ) 


and resubstituting Y!j^ = 1 y^^ + 


iu (i|) ' -2; 


ln(l+d'Z) into eq.uation 4 «5.23j 


it transforms to 


^CiIfz) Cos [01n(l+d«z) ] 


+ D(l+d«Z) Sin [ 01n(l+d»Z) ] - ) Ml+4«Z)}^ (4,3.25) 
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Equation 4»3«25 represents the general solution for equation 4.3,16. 
For -the lower layer the governing differential equation is s 


2 

dZ dZ 


dy!, 

dz" 


'11 


z^izizt 


Suhsti tuting 

(l+diz) = n 

V 

and Ti = e and proceeding in the sane way as hefore equation 4*5*‘I2 
transforms to 5 


+ 2 -J- + (*.-2) _ 

dv dv dv 


11 


+ = 0 .... (4.5.26) 


The general solution of equation 4*5.26 is similar to that 
of equation3,5.23 nnd is given hy (ince 1956> Kamke 1959 aiid Kosko 1965) 


y^l^ = (l+d'Z)“'‘"’ a'Cosh [a»ln(l+d'Z)] Cos[T'ln (l+diz)] 

+ B‘ Cosh[ a'ln(l+d’Z)] Sin [ t' ln(l+d’ z) 3 + C Sinh [ a’ ln(l+d’ z)] 
Cos [t > ln( 1+d’ Z)] +1)’ Sinh [ 0 ’ ln( 1+d' Z) ] Sin[ t' ln(l+d*Z) ] } 

... (4.3.27) 

The arhitrary constants A,B,C,D, A*, B>,C' and D' alongwith 
the unlcnovm load, Q.^-^ or , such that the soil goes into plastic 
region upto a depth of can he found from the four boundary condi- 
tions at the tip and top of the pile and the five oompatihility 
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conditions at the interface. 

For the case (i) i,e. for free-free pile acted upon by axial load 
P and later farce the boundary conditions are ; 

I atxcO i.e, atZtaO 

»«0 . 4 . . (4.5.28 n) 


a) B ,M ss 0 Q 

^ y, 

Bl(x) 

dx" 


u 




or 


u 


0 


dZ 


(4.5.28) 


and 


b) S.P = 

d^Y dy 

dx 


dx 


= *Q.t 


> • *« (4.5^9a) 


5 


dy 


or 


2 

,, -^lu ^ 


1u 


tR ^tD 


.... (4.5.29) 


■where = Uon-Bimensional shear force at the top 

II at X » L.| or Z « the five compa'tibility conditions are i 

a) Deflection oompatibility : 


....(4.5.30) 
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"b) Slope oompatibilitys 

dZ ” dZ 

c) Bending moment oompatibility s 

^ , ilk 

2 . 2 . 
dZ dZ 


d) Shear compatibility s 


5 5 

d V* d yf 

'^11 


dZ' 


dZ 


5 


and 


e) Pressure compatibility 




III at X = L or Z = Z» 


a 


) B.M. = 0 


d yl 


11 


dZ 

and 

b) S.P = 0 


U 


dy’ 


^11 

dZ^ (1+d'Z') dZ^ (l+d'Z’)^ 


4d« ^11 


11 


0 


• • • • (4 *3 *51) 


• • • • (4 *5*32) 


.... (4.3.35) 


.... (4.3.54) 


.... ( 4 . 3 . 35 ) 


.... (4 * 3 • 5^ ) 


[ 
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By using th© aBove .9 iDoundary conditions (equations 4 » 5*27 
to 4«3»55) the constants A,ByC,D, and D’ and the lateral load 

required for a particular value of to go into plastic resistance 
are evaluated. 

Bor the case (ii) i,e. free-free pile acted upon hy a moment 
ilj. and an asdal load P and case (iii) fixed-free pile thq general form 
of the solutions, remain same. The houndary conditions at x = L and 
compatibility conditions at x = also remain same, whereas the 
boundary conditions at the top of the pile, at x = 0 are j 


Case ii) For free-free pile ¥/ith only P and acting at top, 
at X = 0. 


a) B .11 . = 


• . • . (4 *5 •37 


or 


d^yi 


dZ 


IL 


tE 


2 = 


.... (4.3.37) 


where M 
and b) 


tD 


= Eon-dimensional moment 


S.P = 0 


dZ^ 


+ 4d’ 


,2 

tik 

dZ^ 


+ U 



= 0 


.... (4 . 3 * 38 ) 


and for case (iii) Fixed-free end conditions pile vath only P and 
0,^ acting at top 

at X = 0 


a) slope a 0 


i.e, 


dyi 


If 


dZ 


0 . 


f.l. T. 

CENTRAL ItBRARy 


( 4 . 3 . 39 ) 




21 IKS 



52 


b) S .P « - 


cL^y’ 


1f 


dZ 


3 . 


+ 4d' 


az" 


tR ^1f 


where « Non-dimensional shear forse at the top. 


4.4 RESULTS m DISCUSSION. 


.... ( 4 . 5 . 40 ) 


Numerical results are obtained for normal ranges of variables, 

based on the generalized solutions presented in the previous section. 

Here again Z' and c are taken as 4.0 and 0.6 respectively. Results 

are obtained for =0.5 and 0.6 wherein U = 0.91 93 • As most 
cr 

probably only a small portion near the ground surface may go into 
plastic region, the value of Z^ is likely to be around 0.25 to 0.5. 
Hence numerical results are obtained for Z^ = 0,2, O.4 and 0,6 which 
covers a wide range of practical interests. 

Pigs. 4*5 a^ad 4.4 show the effect of Z^ on and 

respectively for rn?- = 0.3 snd 0,6. It is seen from these that, for 
^or 
U 

— ■= 0.3, as Z^ increases from 0.2 to 0.6 the lateral force required 
cr 

or , for a depth of the soil at the top to go into plastic 

region, also increases considerably. Hence as expected in this case 

as the lateral load is increased more and more soil at the top goes 

into plastic zone. Whereas for = 0,6 even for a small 

cr 

increment in the lateral load, the depth upto which the soil at the 
top goes into plastic region considerably increases. This is because 
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at higher axial loads on the pile, the soil at the top does not take 
much lateral load to reach ultimate resistance. 

Fig. 4*5 shows the effect of on the deflection coefficient 

at Z = 0 for 2 — ^ 0,5 and 0.6 . Here it is seen that as the depth 
^or 

of plastic zone increases the deflection coefficient y!j at Z z= 0 

TF 

increases oonsidorahly for both = — =0.3 and 0.6. Fig, 4,6 depicts 

^or 

the effect of Z. on yl at Z=0 for = 0.3 and 0.6 . Here again 

the same trend is seen as in Fig. 4.5* 

Fig, 4*7 shows the relationship between the depth of plastic 
zone, Z^, and the corresponding lateral load required, Q,^^, for fixed 
head pile for - ®*3 snd 0.6. Here again as expected the lateral 

load required at ihe top increases with increasing Z^, 

Fig. 4*8 shov/s the effect of Z^ on fixity moment at the fixed 
end of the pile. It is seen that as Z^ increases the moment at the 
fixed end increases as can be expected. 

Fig, 4*9 shows the e ffect of increasing Z^ on the deflection 
coefficient at top, y’ , for a fixed head pile for u/u = 0»3 end 0,6. 

4 X _ C!l? 

It is seen that the effect of is considerably less on the 

deflection coefficient for fixed head piile. It is also seen virhen 
compared to froG-froG pile the deflection and moment coefficient 
values axe considerably less for fixed-free pile for the same value 
of and lateral load. 



54 


4 .5 CONCLUSIONS 

•As stated earlier the floxaral 'beha:viour of axially and 

laterally loaded piles can be reasonably carried out at working 

loads on the concept of soil modulus. However, v/ith increasing 

lateral loads the soil in the top region reaches its ultimate 

resistance value and as such for any realistic analysis one has to 

take into account the elasto-plastic nature of the soil, With this 

in view the analysis presented in this Chapter can be used for 
/ 

analysing the problem of axially and laterally loaded piles with 
variable cross-section, when the pile is subjected to lateral loads 
which arc greater than worlcLng loads. Also the analysis brings out 
the fact that by assuming the soil to be elasto-plastic and analysing 
the problem, the effect of non-linearity between load-deflection 
curves observed in the field can bo accounted for. 

The results obtained show that the boundary conditions at the 
top depth of plastic zone, and the axial force coining have an over- 
whelming influence on the flexural bohaviour. 

The solutions being obtained in closed form, those can be 
also useful for checking the numerical efficiency of other approxi- 
mate methods which are in general used for solving such complicated 
problems , 





























CHiiPTSR 5 


EROBiffilLISTIC APHOACH TO DESIGN OP PILES. 

5.1 INTRODUCTION 

In the previous chapters it was considered that the soil 
properties', pile properties and the loads coming on the pile are 
detoministic in na'turo. However, in reality they are random in 
na-ture rrith certain control in si%A;xvi sion . The necessity for 
prohahilis-feic considerations in soil mechanics arises from the fact 
that soil is not a material made to satisfy a given specification, 
like steel and concrete, hut a na'tural product of which the properties 
may vary from place to place considerahly. This variation may he 
due to the uncertainties in the soil properties or due to strati- 
fications that contain several distinct layers. The prohahilistic 
analysis of the latter is relatively involved and as such imich 'work 
has not hecn done taking this factor into account. However, many 
authors (hazard 1961, Langejan 1965 » Lumh 1966, V/u and Kraft , 
Resendiz and Herrera 1969, Shuk 1969, Biomatowsky I969, ^‘’^u and 
Kraft 1970 cjnd Benjamin 197 ^) have considered the randomness in the 
soil properties in analysing prohlons of foundation engineering 
and soil mechanics ignoring the spatial correlation. 

Analysis of the random variation of the strength properties 
of materials such as concrete and steel hy many authors 
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(Freudonthal 1947> 1956, I96I jIumlD 1966 and others) show that a 
numher of i^rohability distributions can. he used to describe the 
variation. These include the normal, log-normal and extreme value 
distributions. Since the controlling factors of the distributions 
are not very well known, it is necessary to determine the appropriate 
probability function empirically by fitting to the available 
experimental data. 

Eventhou^ in recent times few studies have been reported 
wherein probabilistic api^roaoh is used for the design of foundations, 
there is paucity of literature regarding the use of this technique 
for the design of pjile foundation. As the design of piles is done on 
understrmgth and overload factor basis (Brons 1964), the concept of 
these factors of safety, from probabilistic approach is considered 
at some length in this chapter. The effect of distribution function, 
’Large sample theory’ and t small sample theory’ (Alder and Eocssler, 
1961) as affecting the factor of safety is brought out. 

5.2 REVISE OP UTEEATBEE 

As stated earlier, there is not roach literature available 
regarding the application of probabilistic approach specifically for 
the design of pile foundations. Eith this in view a, brief review of 
the existing literature on probabilistic approaches to soil mechanios 
problems in general is presented belovr. 


Eventhou^ as early as 1947 Preudenthal enunciated the 
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(Preadon'fchsl 1947 ? 1956 ? I96I ,11111113 1966 and others) show that a. 
nunher of prohahility distrihutions can he used to describe the 
variation# These include the nonaal, log-nonnal and eDctreme value 
distributions# Since the controlling factors of the distributions 
are not very well known, it is necessary to determine the appropriate 
probability function empirically by fitting to the available 
experimental data. 

Evonthou^ in recent times few studies hove been reported 
v/horein probabilistic aiDproooh is used for the design of foundations, 
there is paucity of literature regaining the use of this technique 
for the design of pile foundation. As the design of piles is done on 
understrcaagth and overload factor basis (Brons 1964)? the concept of 
those factors of safety- from pixibabilistio approach is considered 
at somo length in this chapter. The effect of distil bution function, 

* Largo sample theory’ and l small sample theory^ ( ilcler and Rocsslor, 
1961) as affecting the factor of safety is brought out. 

5,2 REVTE?? Of LITEHATUKE 

is stated earlier, there is not rijuch literature available 
regarding the application of probabilistic approach specifically for 
the design of pile foundations. With this in view a brief review of 
the e3dsting literature on prohahilistic approaches to soil mechanics 
problems in general is presented below. 

Bventhough as early as 1947 freudenthal onunoiatod the 
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classical lixcTDabiHstic approach to structural problems, only recently, 
hazard (1961) presented a discussion on probabilistic approach to soil 
ncohanics problem very briefly. 

Lange j an (1965) analysed the problem of the stability of a 
slope \ 7 hcrein .. it is assumed that ihe governing soil property is 
only the shear strength, S^, and is taken to bo a random variable. The 
design shear strength, S^, is taken as equal to where F is the 

factor of safety applied. The probability of failure is stated as 

^f 

^ f Sf t ¥~ ^ wherein = true shear strength of the soil : . 

and = mean shear strength* evaluated from the laboratory and/ 

or field teats. 

Further, it is assumed that the distribution of to be log- 
normal and the probability of failure of the slope for different values 
of factor of safety, P, variance of S^, s, and the number of samples, 

N . It is scon from the results presented that variance and number of 
samples tested have overwhelming influence on the failure probability. 

Lange j on (t965) further studied the economic aspects of the 
failure of an eroavation and arrived at an optimum slope based on the 
statistical-GOonoBliG decision theory. 

Lumb (1966) presented the results of tests conducted on four 
different typos of soils from Hong Kong, Viz.: (a) soft marine clay 
deposited in shallow waters (b) alluvial sandy clay (c) a recidual 
silty sand (d) a recidual clayey silt. Tho properties studied 


I 
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includG Attcrberg limits, groding and for undisturlDed samples, 
strength and compressibility ch£3ractGristics. Herein it is shown that 
the properties follow Gaussian distribution or related to Gaussian 

2 

distribution. Comparison of the observed and theoretical values by 
tost show good aggroement betv/een them. In the analysis only 
■uncertainty in the soil properties is considered* 

Wu and loraft (196?) considered the randomness in both soil 
propcT'tios and loads in the analysis of foundation safety. The experi- 
mental data used in the analysis v^as obtained from a total of about 
twenty five borings in an area of about 15 OOO sq..ft. and about 100 
unconfinod •bests wore made on 2-inch diameter "bube samples. The site 
consists of a brown sandy clcy 'till and a, grey clay till with lenses 
of fine to medium sand. The results of the above "bests fit the log- 
normal distribution fairly 'well. The results of ihe. standard penotra- 
tion tests conducted on sands in 12 borings also fi'b the log-normal 
dis’bribution reasonably good, 

Y/u and Craft (19^7) further verified Tersaghi and Peck’s 
( 1948 ) recommendation which says that, the noasurod strength with in a 
significant depth from each boring be averaged and the smallest 
average bo used for design purposes. It is seen that if the minijaun 
average strength is used for the design then the variations in soil 
properties and the size of failure area do not exert a predominant 
influence on safety. In the above oalculations,tlie dead load is 
assumed to bo dotorministic and the live load is assumed to be 
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normally distri’butod, and the wind load to be eDqDonentially distributed. 
It is also shown that an increase in number of borings in a given site 
improves the probability of safety considerably whereas an increase in 
number of samples in each boring does not have much effect on the 
probability of safety, 

Ifiihn and Burton (196?) proposed a statistical method to control 
the "dry density for a rood construction." In the proposed sampling 
scheme a method is outlined whereby engineering decisions can bo taken 
on the number of somplos to be tested to satisfy various engineering 
roquiremonts such as the error in the tost results and the desired 
limits of accuracy for a specified probability. It is seem from the 
tost results presented that the dry density fits the normal distribution 
fairly well. 

Rcsondiz and Herrera (1969) studied the sottlemoit Eaa 4 rotation 
of rectangular foundations on randomly compressible soils for both 
flexible and rigid foundations . Herein all variations in compreasi- 
bility ocouring in thoh rizontal direction are random. It is shown 
that the average settlement and rotation of a foundation can he 
regarded as normally distributed also other parameters being equal, 
the probability of rotation exceeding a certain value increases with 
the rigidity of the foundation. 

Shuk (1969) presented an analysis of variability of natural 
soils in the design of footings with regard to settlements due to 



60 


oonsolidaljion. The "total ejqpectod luildiug cost) "which is the suia of 
initial cost plus "the ejjpeoted maintain once cost^ for the conditional 
prolaalaility of exceeding the gi-von maxinum allowable total or differen - 
tia.1 settlement is plotted against the probcvbility of exceeding a 
maxiinuja allowable settlement. The probability corresponding to the 
miniimm total expected cost gives the "Optiinam design probability". 

BiematoWskL (1969) made an analysis of stability of slopes 
•with probabilistic considerations. In the analysis all the parameters are 
taken to be random variables except the geometrical dimensions. 

Analysis is prcaccated for both normal and log-normal dis"fcribution of 
random variables and for different degrees of slope. 

Benjamin (1970) bhe first time used Bayasicn docieion tlio-jry 
approach in predicting the probable settlements of foundations. Herein 
an apriori dis'tribution is assumed before the collection of data bared 
on the oxporionco and professional judgement. After filling the data, 
Pootcrior distaxLbu"tion is obtained by using Bayesian theory* The total 
expected loss XThich is the sum of sampling cost and expected terminal 
loss, is minimized "to get the optimal decision. 

Wu and Kraft ('V970) presented the analysis for the probability 
of failure of slopes. For various values of conventional factor of 
safety the failure probability has been evaluated and statistical 
decision theory is applied to the selection of optiraam slope. Here 
again the optinum design is one v/hich minimizes the expected cost or 
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nuliiiDiizGs ijhG ca^ectod utility with in. the Gcoxiomic or physical 
constraints imposed hy Ihe problem. It is shown in this work that 
the optimum safety factor is not unique quantity hut depends on 
the stability number of Ihe slope, 

Kogan and Lupasko (l970) presented the variational method 
of evaluating slope stability at limiting equilibrium conditions. 

Here a general shape of sliding surface is assumed, whereas 
eorlica: works were confined to circular arcs, logarithmic spirals, 
etc. 

Ignatov ( 1970 ) used the method of end differences to solve 
the problem of designing rectangular slabs resting on a statistically 
non-homogeneous base and subjected to uniformly distributed loading 
applied over rootangular areas. 

5,5 M/JiYSIS 

As stated earlier in general the design of laterally loaded 
piles is carried on the basis of over-load and under- strength factors 
(Broms, 1964a, 1964 • Thus the failure of Ihe pile can take place 

either' when (i) the actual loads coming onto the structuro exceed 
considerably than those used while designing and/or (ii) if the 
strength paxametors of the soil has been over estimated. Soil being 
a natural product has considerable variation and as such one should 
really think in terms of probabilistic variation of strength while 
dealing with the design of laterally loaded piles at ultimate 
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resistance. Having recognized the fact that the strength should 
ho taken as a random function with n certain prohabilistlc distrihution 
based on the experimental data, one can then consider the partial 
factor of safety in terms of uncertainties. In ishat follows, it is 
assumed that all other factors arc doteiministic whereas only the 
properties of tho soil ore considered with uncertainty. However, for 
a complete analysis one has to consider the uncertainties both in 
loadings as well as in properties of the soil and pile material. Also, 
it should ho noted that the following discussion deals only with the 
factor of Swafetjr-concept and does not give complete analysis of the 
problem of design of piles for resisting la,teral loads. 

For a short, rigid pile in a cohesive soil, it is assumed that 
the failure talces place, when all along the embedded length of the pile— <• 5^1 
the soil reaches tho ultimate resistonoe which is equal to where 

Cl, a^nconfined compressive strength, and d = diameter of the pile, 

^ a 

Lot tho shear strength be tho random voriable, with 
certain distribution. Then in the conventional method, the design 
shear strength, is given by s 

, S/F 

in which «= mean value of the shear strength obtained from 

laboratory and/or field tests, 
and F » Conventional factor of safety. 
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How talcing uncertainties in soil strength, the prohahility 
of soil reaching ultimate strength can he expressed as follows t 

Pf - pr{Sft- ) < 0] .... (5.1) 

in which = the prohahility of soil reaching ultimate resistance 
Q 2 id S|,.j. =» True in - aitu shctix ctrength 

If the number of samples tested from a population axe loss 
than jOy then small sample theory, given hy Alder and Eoessler 
has to he used, wherein Student^ s jj-distrihution is used to 
approximate the actual distribution 'of S^, 

Case (i) If is normally distributed, then 



/l(Sf • Sf)^ 

in which s e ^ =;r-r * whore IT = number of saonplos tested. 

iwl 

t « the variable of Student's distribution function. 

Substituting for S^.(. from equation 5.2 into equation 5.1 , it reduces 
to : _ 

Pf - P>[t > ~ ^ 1 (5.5) 

In a given case the values of S^., s, P and IT are knov/n ; then P^ 
can he obtained from tables of the t-distribution found in many 
statistical handbooks (Pisher and Yatesr 



TiiBIiHl I . TalUGs of for » 0,8 Kg3,/cia‘^ and for different valacs of Number of 

Tests (N), Standard Deviation (s), and Pan tor of Safety (P) t 



6d 
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TalDlo I gives the magnitude of calculated from 

values of IT, c, and P that represent the range of practical possibilities. 
Hoi^ein is assumed- to be 0,8 Kgs./sq..om. Prom this table it 

appears, that if a pile is designed with a given value of the safety 
factor P, the probability of soil reaching ultimate resistance may 
vary within wide range of values, Also it is seen that there is a sharp 
increase in for small values of IT. 

Case (2) If is log-normally distributed s 

Now the = ,P [ • |i~- <0 J can be written as 

p£ » P P< Sj3 .... (54) 

Since the quantitios on either side are not nogail«e, talcing logarithm 
on both sides it tr«ansforms to : 

= . P [ la(S^^) < ln( S^) ] .... (5 .5) 

It is well known that, G. M •••• (5*6) 

in which JLM <■ arithmatio mean of a set of numhers. 

G.M = Geometric mean of a set of numbers. 

If possible values of then from 

equation (5 ♦6) 

l/!t 

Csj, + Sf2 + ..-.j- Sgj) ^ w .... (5.7) 
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Equality holds good only if e 


1 .e. 


- 1 ^ ■ 

1 ^ to ( iLi'^fi^ 

to(S^) 1 toXsp' , 


» » • • (5 *6 ') 

....(5.9) 


From equation 5*5 equation 5*9 it is seen that only a lower hound 
can he given : 

- P [to(S^^ P) < ln(sp 1 


P[ ln(s^/) 


»*.. (5.10) 


Again hy siaall aanple theory 


t to - In 




.... 


in which 


/i7i 


_ S (to Sj - In S^^)‘ 


Now suhstituting for ln(Sj^^) from equation 5-11 into equation 5*10 
la ohtednod as : 


P^> P'it 


> " In P ] 

a 


.... (5.12) 


For different values of N, s and F ,thG lower hound for 


is evaluated. 
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TaTsle II repujGseiats values for different values of 
11,3 and P. Here a^odn it is seen that for a given factor of safety 
the prohability of failure varies widely. It is also seen from the 
Tables I and II that the is considerably less for log-nonaal 

distribution of when compared to that of normal distribution of 
for the same values of IT, s and P. 

Usually a sample from a population is said to be ' Large < 
if the mmber of samples tested is greater than 50 * However it should be 
noted that there is no clear cut boundary between ’Large Sample ’ and 
' Small Sample* . 

According to central limit theorem "if X has a distri - 
bution with moan, m,. and Standard deviation, ® , and if all possible 
samples of H variates axe drawn then the sample mean, X, will have 
a distribution which for larger and larger N, approaches the normal 
distribution with mean, m and Standard deviation a / ” (Alder and 

Rocssler, 1961 ) . 

i.e. Z » ^ N(0,1) for large N *...(5.13) 

‘^X 

where ^ m o / ^ , 

So it is seen from this theorem that whatever may be the 
distribution of the sample mean S^ will be normally distributed 

for large N. V/lth this in view a general equation is given below for 

Ip • 
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^ [Sft < / 3 (5.14) 

■but according to the above theorem 

Z - N( 0 ,l} for large N. ( 5 .I 5 ) 

s 

sZ 

%t “ + ^f •.*• (5.46) 

Substitutii^g ior ili equation 5*14 -ixom equation it is seen, 

that, 

Ptz>^ 3 ....(5.17) 

Using equation 5«47 ^”ox different values of U, s and P , P^ is 

O 

caloulated for s 0.8 Kgs/cm*^ and is presented in Table III* 

Here again it is seen that the effect of N,s end P is quite 


I # 


considerable . 



Ti’vBlE III. Vaiues of for is 0.8 Kgs/cm and for different values of Karaber of tests (1T), 
Standccrd Deviation (s ) and Factor of Safety (P) . 
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5.4 Suac-ESTIOUS EOR further EESEiUlCH 

Baaed on the review of literature and the above 
analysis certain suggestions for further research are made , 

The above analysis clearly brings out the fact that in the actual 
design it is worthwhile to consider the probabilistic variation of 
soil properties which would give the probability of failure rather 
than using arbitrary factors of safety. Also the analysis shows 
the effect of distribution functions, large and small sample theory 
as affecting the probabilistic analysis. 

However for a realistic analysis, the loads and pile 
properties are also to be considered as probabilistic rather than 
deterministic. Also regression analysis is to be used to consider 
the spatial correlation. The actual distribution function of the 
soil property should be found out by fitting empirically to the 
experimental data available. In the absence of data subjective 
probabilities have to be used in making decisions. With all this 
in view, it is seen that, there is ample scope for further research 
in probabilistic approach to soil mechanics problems. 
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